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Abstract 

We consider a linear Hamiltonian system consisting of an infinite heat reservoir and a 
classical particle. The initial data are supposed to be a random function which has some 
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1 Introduction 



The paper concerns problems of long-time convergence to an equilibrium distribution for a 
coupled system consisting of a field and a particle. For one- dimensional chains of harmonic 
oscillators, the results have been established by Spohn and Lebowitz in [28] , and by Boldrighini 
and others in [2j. Ergodic properties of one-dimensional chains of anharmonic oscillators cou- 
pled to heat baths were studied by Jaksic, Fillet and others (see, e.g., |211[T1]). In [6| W\ [8| [TO], 
we began to study the convergence to equilibrium for the systems described by partial differ- 
ential equations. Later on, the similar results were obtained in [9] for ci-dimensional harmonic 
crystals with d > 1, and in [IT] for a scalar field coupled to a harmonic crystal. 

Here we treat the linear Hamiltonian system consisting of the scalar wave or Klein-Gordon 
field (p{x), X G W^, coupled to a classical particle with position in g G M'^. The Hamiltonian 
functional reads 

Hi(p,ir,q,p) = HAiq,p) + Hb^V^t^) + <? " {^V,p), (1-1) 

where 

HA{q.p) = \[\p\' + u'\q\'), 

HBi^,n) = ^l{\V^ix)\' + m'\^ix)\' + \nix)\')dx. 

Here m > 0, u > 0, (p{x) and tt{x) are real-valued functions, "■" stands for the standard 
Euclidean scalar product in M'^, (-, ■) denotes the inner product in the real Hilbert space L'^{W^) 
(or its extensions), q,p eM.'^, d > 3, and d is odd if m = 0. 

We assume that the initial data Yq := (v^o? ttq, go^Po) are a random element of a real func- 
tional space S consisting of states with a finite local energy, see Definition 12.11 below. The 
distribution of Yq is a probability measure /iq of mean zero satisfying conditions S1-S3 below. 
In particular, we assume that the initial measure /zq satisfies a mixing condition. Roughly 
speaking, it means that 

Yo{x) and Yq^i/) are asymptotically independent as \x — y\ ^ oo. 

We study the distributions /if, t G M, of the random solution Yt := {ipt,''^t,qt,Pt) at time 
moments t G M. Our main objective is to prove the weak convergence of the measures fit to 
an equilibrium measure /Zoo, 

/it ^ Atoo, t -> oo, (1.2) 

where the limit measure //oo is a translation- invariant Gaussian measure on S. We derive the 
explicit formulas for the limiting correlation functions of /Zoo- The similar convergence holds 
for t — )■ — oo since our system is time-reversible. 

Let us outline the strategy of the proof. When the bath variables {ipt, vr^) are eliminated 
from the equations of the coupled system, the particle evolves according to a linear Volterra 
integro-differential equation of a form (see Eqn (13. 2p ) 

t 

qt = -co^qt + j ds + F{t), t G M, (1.3) 
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where D(t) is a matrix- valued function depending on p, F(t) is a vector- valued function de- 
pending on (990, ttq). Therefore, our first objective is to study the long-time behavior of the 
solutions to Eqn.f ll.3]) (see Theorem 13.11 and Lemma [3.61 below). The deterministic dynamics 
of Eqn (II. 3p has been extensively investigated by many authors under some restrictions on the 
functions F{t) and D{t): Myshkis [21], Driver |1] and others (for more details on the results 
and problems in the theory of equations with delay, see, e.g., the survey paper by Corduneanu 
and Lakshmikantham [3J, the monograph by Gripenberg, Londen and Staffans ^16j and the 
references therein). 

The linear stochastic Volterra equations of convolution type have been treated also by 
many authors, see, e.g., the papers by Appleby and Reynolds [1], by Karczewska [22] and the 
references therein. In the papers mentioned above, the solutions of Eqn (II. 3p were studied 
assuming that F{t) is a Gaussian with noise and the memory kernel D{t) has the exponential 
decay. However, in our case, F{t) is not Gaussian white-noise, in general. Moreover, if m 7^ 0, 
the decay of D{t) is as (1 + t)"*^/^ (see Appendix B). 

In recent years the nonlinear generalized Langevin equation, i.e., the equation of a form (cf 
Eqn.dmSD) 

t 

qt = -"^Viqt) - y r(t - s)qs ds + F{t), t G M, (1.4) 


with a stationary Gaussian process F{t) and with a smooth (confining or periodic) potential 
V{q), has been investigated also extensively, see, e.g., [201 ISSl [271 EO]- In particular, the 
ergodic properties of (ll.4p were studied by Jaksic and Fillet in [20], the qualitative properties 
of solutions to Eqn. (ll.4| were established by Ottobre and Pavliotis in [25]. 

In this paper, we study a linear "field-particle" model. However, we do not assume that 
the initial distribution of the system is a Gibbs measure either absolutely continuous with 
respect to a Gibbs measure or "close to equilibrium". Therefore, in particular, the force F{t) 
in Eqn. (ll.3| is non-Gaussian, in general. 

The key step in our proof is the derivation of the asymptotic behavior of the solutions 
Yt for the coupled field-particle system. We prove the following asymptotics in mean (see 
Proposition 15. ip : 

{Yt,Z)r^{Wtiipo,'iTo),niZ)), t->oo, (1.5) 

where Wt is a solving operator to the Cauchy problem for the free wave or Klein-Gordon 
equation (I2.10p . and the function Q{Z) is defined in (I2.17p . This asymptotics allows us to 
apply the technique from [HI [10] , where the weak convergence of the statistical solutions has 
been proved for free wave and Klein-Gordon equations. Namely, we divide the proof of (II. 2p 
into two steps: we first establish the weak compactness of the measures family {pt, t G M} (see 
Section [1|), and then we prove the convergence of the characteristic functionals of the measures 
fit (Section [6]). 

In conclusion, note that the similar results remain true for a linear Hamiltonian system 
consisting of heat reservoirs coupled to a single particle. In this case, Hamiltonian reads 

N N 

^HB{^i,'ni) + HA{q,p) + q-^{V^k,pk)- 

i=l k=l 

Also, convergence (11.21) holds for a system consisting of a field if{x), x G M'^, and an oscillator 
with position g G M\ if Hamiltonian of the system equals Hb{(P, vt) + HA{q,p) + qCVip, Vp). 
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The paper is organized as follows. In Section 2 we describe the model, impose the conditions 
on the coupled function p and on the initial measures /io and state the main results. Section [3] 
is devoted to the limit behavior of solutions of Eqn fll.Sp . In SectionlHwe prove the compactness 
of the measures family {/i*}. The asymptotics (II. 5p is proved in Section [5l In Section |6] we 
complete the proof of the main result and establish the convergence of characteristic functionals 
of Ht to a limit. In Section [7] we study the ergodic properties of the limit measures /loo- 
In Appendix A we prove the existence and uniqueness of solutions of the coupled system. 
Appendix B is devoted to the asymptotic behavior of the solutions of (II. 3p in the case when 
m 7^ 0. Appendix C concerns the case of non translation-invariant initial measures. In 
Appendix D we give an application of the results to the case of Gibbs initial measures. 

2 Main Results 
2.1 Model 

For simplicity of exposition, we consider the case c? = 3 only. Taking formally variational 
derivatives in (II. ip . we derive the following equations of motion, 

0t(x) = 7rt(x), 7i-t(x) = (A-m2)(^j(x) + gi- Vp(x), x G M^ t G 
(It = Pt, Pt = -(^"^Qt + j Vp{x)ipt{x) dx. 



(2.1) 



For the system (12. ip . we study the Cauchy problem with the initial data 

Vt{x)\t=Q = (po{x), 7it{x)\t=o = vro(x), X e M^, qt\t=o = go, Pt\t=o = Po- (2.2) 

Let us write (pt = {<^t{-),TTt{-)), 6 = {(lt,Pt), Yt = {4>t,^t)- Then the system ([2lT])-([22]) writes 
as 

Yt = C{Yt), teR; Yt\t=o = Yo. (2.3) 

Now we formulate the conditions on p, u and m. Let Km = (-^m.ij )i*j=i stand for the 3x3 
matrix with the matrix elements Km,ij, 

Km,,-=i2.)-^l^§f^dk, m>0. (2.4) 

Al. The matrix (cu^ — m^)J — Kq is positive definite, where I denotes the identity matrix. 
Al'. The matrix cj^J — is positive definite. 

Note that if m 7^ 0, then condition Al implies Al'. Otherwise, if m = 0, conditions Al 
and Al' coincide. 

A2. The function p{x) is a real- valued smooth function, p{—x) = p{x), p{x) = for |a;| > Rp. 
A3. For any keR^\ {0}, p{k) = J e^'^Xx) dx 0. 

2.2 Phase space for the coupled system 

Let II ■ II denote the norm in L^(M^). We introduce a phase space S. 
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Definition 2.1 (i) Let us choose a function ({x) G C^(M^) with ({0) 7^ 0. Denote by 
H^^^iM.^), s G M, the local Sobolev spaces, i.e., the Frechet spaces of distributions ip G D'(R^) 
with the finite seminorms \\^\\s,r '■= ||A^ {({x/R)ip) \\, where A'^ stands for the pseudodifferential 
operator with the symbol (k)^ , i.e., 

A'^ := F,-U{krm), (k) := vW+l, 
and if) is the Fourier transform of a tempered distribution ip. 

(a) Ti = Hl^^{M.^) © i7j°^(R'^) is the Frechet space of pairs = 7r(x)) with real valued 

functions (f{x) and tt{x), which is endowed with the local energy seminorms 




(|(^(a;)p + |V<^(x)|2 + \7i{x)\'^)dx < 00, Vi? > 0. 



|x|<i? 

(Hi) S = "H © © is the Frechet space of vectors Y = {(f){x),q,p), with the local energy 
seminorms 

\\Y\\ln=ml+\q\' + \p\', Vi?>0. (2.5) 
(iv) Fors G R, denote W = ///„|'(M^) © i/f„,(R^), S' = H'®M.^®M.^. In particular, K = H°, 

Using standard techniques of pseudodifferential operators and Sobolev's Theorem (see, e.g., 
[T7]). it is possible to prove that £^ = £ C £~^ for every £ > 0, and the embedding is compact. 



Proposition 2.2 Let conditions Al' and A2 hold. Then (i) for every Yq G £, the Cauchy 
problem Ii2.3\) has a unique solution Yt G C{W,£). 

(a) For every t G M, the operator St : Yq Yt is continuous on £. Moreover, for every 
R> Rp,T>0, 

sup \\StYo\\e,R < C{T)\\Yo\\£,R+T. (2.6) 

\t\<T 

Proposition 12.21 is proved in Appendix A. 

2.3 Conditions on the initial measure 

Let {Q, S, P) be a probability space with expectation E and B{£) denote the Borel a-algebra 
in £. We assume that Yq = Yq{u,x) in (12. 3p is a measurable random function with values in 
{£, B{£)). In other words, (w, x) 1— )■ Yq^uj, x) is a measurable map x — )■ with respect to 
the (completed) a-algebra S x i3(]R^) and i3(]R^). Then Yt = StY^ is also a measurable random 
function with values in {£,B{£)) owing to Proposition 12.21 We denote by /io('^^o) a Borel 
probability measure in £ giving the distribution of Yq. Without loss of generahty, we assume 
(f2, S,P) = {£,B{£), fio) and Yq{ijj,x) = u{x) for fio^du) x dx-aAmost all {ijJ,x) G ^ x R^. 

Set V = Vo®M.^® M^ Vo := [C^{R'^)]\ and 

(Y, Z) := {(p, f) + q ■ u + p ■ V for F = (0, g,p) G and Z = {f,u,v) E V. 
For a probability measure fi on £ denote by /t the characteristic functional (Fourier transform) 

il{Z) = [ exp(i(r, Z)) ix{dY), ZeV. 
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A measure /z is called Gaussian (with zero expectation) if its characteristic functional has the 

form fi{Z) = exp{ — -Q{Z,Z)}, Z G P, where Q is a real nonnegative quadratic form in P. 

A measure fi is called translation-invariant if fi{ThB) = fi{B), \/B G S{S), h G M^, where 
ThY{x) = Y{x - h). 

We assume that the initial measure /iq has the following properties S0-S3. 

50 /io has zero expectation value, WfQ{x) = J Yq{x) fiQ{dYo) = 0, x E M^. 

51 fiQ has finite mean energy density, i.e., Ed^oP + boP) < oo, and 

E(\ifo{x)\^ + \VMx)\^ + \Mx)\^) < eo < oo. (2.7) 

Denote by fj,Q := PfiQ, where P : {4>o,Qo,Po) G £^ — ?■ 0o G "H. Now we impose condi- 
tions S2 and S3 on the measure . For simplicity of exposition, we assume that (Iq has 
translation-invariant correlation matrices (non translation-invariant case of fiQ is considered in 
Appendix C). 

52 The correlation functions of the measure /i^, 

Qi^{x, y) ■= j (0^(x) ® (t)i{y)) fi^id^o), x,y e M^ 0o = (0o, 0o) = iVo, vto), 

are translation-invariant, i.e., QQ{x,y) = {x — y), i,j = 0, 1. 
Now we formulate a mixing condition for the measure fi^. 

Let 0{r) denote the set of all pairs of open convex subsets A,BcM.^ of distance d{A, B) > 
r, and let cr{A) be the a-algebra in H generated by the linear functional (j) (0, /), where 
/ G [C^(M^)]^ with supp / C A. Define the Ibragimov mixing coefficient of a probability 
measure fi^ on "H by the rule (cf [HI Def. 17.2.2]) 

,(r)^ sup sup l.f(^nB)-.f(W)l ,2.8) 

(AB)60(r) ^ea(^),SGc7(B) ^0(5) 
/iO^(i3) > 



Definition 2.3 VKe sa?/ that the measure fiQ satisfies the strong uniform Ibragimov mixing 
condition if ip{r) — > as r ^ oo. 

S3 The measure yU^ satisfies the strong uniform Ibragimov mixing condition, and 

+ 00 

j r^ip^/^{r)dr < oo. (2.9) 





Remark 2.4 (i) The examples of the measures fiQ with zero mean satisfying conditions 02. 7p . 
S2 and S3 are given in [6]. 

(ii) Instead of the strong uniform Ibragimov mixing condition, it suffices to assume the 
uniform Rosenblatt mixing condition [26] together with a higher degree (> 2) in the bound 
(12. 7p . i.e., to assume that there exists a. 6, 6 > 0, such that 

E (|^o(x)r+' + \VM^)\'^' + ko(x)P+') < oo. 
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r+oo 

In this case, the condition f l2.9p needs the following modification: / r'^a^{r)dr < oo, where 

Jo 

p = mm{6/ (2 + 5), 1/2), a{r) is the Rosenblatt mixing coefficient defined as in (I2.8P but without 
IJ^q{B) in the denominator. 

2.4 Convergence to equilibrium for heat reservoir 

We first study the Cauchy problem for the free wave (or Klein-Gordon) equation, 

r ^t{x) = (A - m'')^t{x), t G R, X G M^, ,^ 
\ ipt{x)\t=Q = ^o{x), 0f(x)|t=o = 7ro(x), 

where m > 0. The following lemma is proved in [6l p. 7] for m = and in [3, p. 1225] for m 7^ 0. 

Lemma 2.5 (i) For any 0o = (v^o? ttq) ^ "H, t/iere exists a unique solution (pt = {(ft{x), 'ft{x)) G 
C{R,n) to the Cauchy problem [KTI^) . 

(a) For any t G M, the operator Wt : (po ^ (pt is continuous in "H. 
Remark 2.6 Denoting (pt = (v^t, T^t), t E M, we rewrite /12.1(J\) in the form 

^t = CB{(pt), teR, <pt\t=o = <Po, (2.11) 

with = ^ ^ 2 0^''^'^ Fourier transform representation, system l{2.10\) becomes 
(pt{k) = CB{k)(pt{k), hence (pt{k) = Qt{k)(pQ{k), where Qt{k) = exp(CBik)t) . Here we denote 

n T \ / sinwif: 

^B{k) = { , , \ , m = ^ 1 , (2.12) 



-u?{k) 



-bj sin bjt cos bjt 



where u = uj{k) = ^/\k^~+^v?. Hence, the solution of l{2.11\) is (pt = Wt(po = Gt{-) * (po, where 
the Green function Qt{x) = Ff7^^[Qt{k)], * stands for the convolution. 

Let Sm{x) be the fundamental solution of the operator —A + m^, i.e., (—A + m?)£m{x) = 
d{x). Since d = 3, £m{x) = ^-^^^ Define, for almost all x, ?/ G M'^, the matrix-valued function 

Q^ix,y) = Qoc^ix-y), where 

- (gooj„=o - 2 I _ ^01 ^11 + + ^2)^00 I • l^-i^J 



Here q^, i,j = 0, 1, are correlation functions of /i^ (see condition S2), and * stands for the 
convolution of distributions. We can rewrite goo in the Fourier transform as 

qooik) = \[qoik)+Cik)qoik)C^ik)y (2.14) 

where ( )^ denotes a matrix transposition, and 
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Remark 2.7 (i) According to [THl Lemma 17.2.3], the derivatives D^q^ are bounded by mixing 
coefficient: Va G 1?^ with |a| < 2 — i — j (including a = 0), i,j = 0, 1, 

\D%'{z)\<Cip^/Wz\), yzeR\ 

Therefore, D'^q^^ G Lp{R^), p > 1 (see ^ p. 16]). Hence, G L\R^) if m ^ by (l27[3ll . If 

m = 0, (12.91) imphes the existence of the convolution £ * q^^. 

(ii) The matrix qoo{k) = (q'J^if^))] j=o satisfies the equilibrium condition, i.e., 

C{k) = co'{k)€ik), C(fc) = -tik). Moreover, g^O(fc) > 0, q'^ik)* = -q'^J^ik). 

Let S{R^) be the set of all functions / G (:7°^(R3) such that sup^ \{x)l^D''f{x)\ < oo for all 
indices a, (3. Denote by Q^{f, /) the real quadratic form on Sq = [S'(M^)]^ defined by 

QZif, f) = {QZi^, y), fi^) ® f{y)) = {Qoo{x - y), f{x) ® f{y)). (2.16) 



Theorem 2.8 (see JE, \^) Let the measure (Iq have zero mean and satisfy conditions {2.1), 
S2 and S3. Then (i) the measures /if = Hq weakly converge as t ^ oo on the space Ti'^ 
for each e > 0. This means the convergence 



/(0)/if(#) J f\<P)l^lo{d(p) as t^oo 
for any bounded continuous functional f{(f)) on "H^^. 

(ii) The limit measure fi^ is a translation-invariant Gaussian measure on "H. 

(Hi) The characteristic functional of ji^ is of the form ft^{f) = exp { — ^Qf,(/, /)}, f E Sq. 

(iv) The limit measure /i^ is invariant, i.e., W^p^ = fi^, t G M. 

(v) The correlation matrices of fif converge to a limit, i.e., for any f E Sq, 

I |(0,/)|Vf (#) = I |(m0,/)|V?(#) ^ QZifJ) a. t -> oo. 

(vi) The flow Wt is mixing w.r.t. fi^ (see also Section^. 

2.5 Convergence to equilibrium for the coupled system 

To formulate the main result for the coupled system we introduce the following notations. Let 
Wl denote the operator adjoint to Wt: 

(0, Wif) = {Wt^, /), for / G [S(R')]\ ^en, teR. 

For Z = {f,u,v) G V, i.e., / G [^^^(M^^)]^, {u,v) G x we define 

n{Z) := f,{x) + a{x) ■ u + (5{x) ■ v. (2.17) 

Here a = (ai, aa, 03),/? = (/3i, /32, /^s) = ^jOt, 

ai{x) := ^ Mir{s)W'_,Vrpods, with po := (p, 0), (2.18) 

r=l { 
3 

(3iix) := ^ j Uir{s)W'_,Vrp'ds, with p° := (0,p), z = 1,2,3, (2.19) 







the matrix-valued function Af{s) = (A/'ir.(s))f s > 0, is the inverse Laplace transform of 
A^(A) defined in ( KWf . and 



3 +~ 



Mx) := fix) + J2 [ (x) {WsV^p\ f) ds. (2.20) 

i=i { 

Definition 2.9 pt is a Borel probability measure in £ which gives the distribution of Yf: 
fitiB) = fioiSi'B), \/B e B{8), t e R. 

Our main result is as follows. 

Theorem 2.10 Let conditions A1-A3 and S0-S3 hold. Then 

(i) the measures fit weakly converge in the Frechet spaces for each e > 0, 

fit — 7 /Uoo as t — )■ oo, (2.21) 
where fioo is a limit measure on S. This means the convergence 

j f{Y)fit{dY) ^ j f{Y)fi^{dY) as t ^ oo 
for any bounded continuous functional f{Y) on S^^ . 

(a) The limit measure fi^o is a Gaussian equilibrium measure on S. The limit characteristic 
functional has the form fioo{Z) = exp{ — -Qoo{Z,Z)}, Z E V. Qoo{Z,Z) denotes the real 
quadratic form on V, 

Q^iZ,Z) = QZmZ),niZ)), (2.22) 

where is defined in Ii2.16\) . and Q{Z) is defined in \2.11^ . 

(Hi) The correlation functions of fit converge to a limit, i.e., for any Z eT), 

j \{Y,Z)\^fXt{dY) =mXt,Z)\^ ^ Q^{Z,Z) ast^oo. 

(iv) The measure fXoo is invariant, i.e., Stfioo = fJ'oo, t G M. 

(v) The flow St is mixing w.r.t. fi^o, i.e., 'if,g G L^(£^,/ioo) the following convergence holds, 
lim / f{StY)g{Y)fji^{dY) = [ f{Y)fji^{dY) ! g{Y)fi^{dY). 



t—^oo 



The assertions (i), (ii) of Theorem 12. lOl follow from Propositions 12.11] and [2712] below. 

Proposition 2.11 The family of the measures {fit, t > 0}, is weakly compact in £~'^ with any 
e > 0. 

Proposition 2.12 For any Z eD , 

fLt{Z) = J exp{t{Y,Z))fit{dY)^exp{-^Q^{Z,Z)}, t ^ cx). 

Proposition 12.111 (Proposition I2.12p provides the existence (resp. uniqueness) of the limit 
measure fioo- Proposition 12. Ill is proved in Section HJ Proposition 12. 121 and the assertion (iii) of 
Theorem 12.101 are proved in Section O Theorem 12.101 (iv) follows from fl2.2ip since the group 
St is continuous in S by Proposition 12.21 (ii). The assertion (v) is proved in Section 7. 
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3 Long-time behavior of the solutions 

Using the operator Wt, we rewrite the system ( 12. ip in the form 



t 



4>t = Wt(t)o + j qs-Wt-sVpUs, xeM^ teM, (3.1) 



t 

qt = -Lo'^qt + (Vpo, ^t) = -oo^qt + j D{t - s)q, ds + F{t), (3.2) 



where 0f = {ipt{-) , 7it{-)) , F(t) denotes the vector-valued function, F(t) = {Vpo,Wt(po), D{t) 
stands for the matrix-valued function with entries 

A,(t) := (V,po, W^iV,/) = {W,p,gt * V,p) = (27r)-3 / hk^ gt{k)\p{k)\^ dk, (3.3) 

Jr3 



with gtix) := F^X^ 



sina;(A;)t 



X G M^ uj{k) = ^/\kfT^, i,j = 1, 2, 3. 



In sections [3] and [5], we study the long-time behavior of the solutions Yt = {(j)t, ^t) of problem 
(12. 3p by the following way. In Section 13.11 we prove the time decay for the solutions qt of (13. 2p 
with F = (see Theorem 13. ip . Then we establish the time decay for the solutions Yt of (12. 3p 
in the case when the initial data (poix) = for \x\ > Rq (Section 13. 2p . Finally, for any initial 
data Yq G £, we derive the long-time asymptotics of the solution Yt in the mean (Section 5). 
At first, consider the Cauchy problem for Eqn. (13.20 with F(t) = 0, i.e.. 



t 

qt = -Lo\t + j ^)^- ds, t G M, (3.4) 



qt\t=o = qo, qt\t=o=Po- (3.5) 
For the solutions qt of (I3.4p -( l3l5|) . the following assertion holds. 

Theorem 3.1 Let conditions A1-A3 hold. Then qt G C^(R), and the following bound holds, 
\qt\ + \qt\ < Cem{t){\qo\ + \po\), where 

, . _ j e""^!*! with a S > 0, if m = 0, , . 

In next subsection, we prove the exponential decay for qt in the case m = 0. The case 
m 7^ is considered in Appendix B. 

3.1 Time decay of qf. Case m = and F = 

To prove Theorem l3.1l we solve the Cauchy problem (I3.4p - (l3.5p by using the Laplace transform, 

+00 



g(A) = j e'^^qt dt, Re A > 0. 



Then Eqn. (13.40 becomes 

X^q{X) = -uj^q{\) + ^(A)g(A) + po + Ago, (3.7) 
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where the entries of D{X) are 

D^j{X)= [ V,p{y){~A + m' + X')-'V,p{y)dy, ^,j = 1,2,3. (3.8) 

We rewrite Eqn. (13. 7p as 

g(A) = [(A^ + uj')I - DiX)] ' (po + Ago) = AfiX)iPo + Mo), (3.9) 

where A/'(A) stands for the 3x3 matrix of the form 

AfiX) = A-\X), with A{X) := {X^ + uj^) I - D{X) for Re A >0. (3.10) 
We first estabhsh properties of A{X). Write := {A G C : Re A > /3} for /3 G M. 

Lemma 3.2 Let m = 0, and let conditions A1-A3 hold. Then 

(i) A{X) admits an analytic continuation to C; 

(ii) for every /3 > > such that v ■ A{X)v > C|A|>p for X G C.p with \X\ > and 
every f G . 

(Hi) There exists a 5 > such that v ■ A{X)v ^ for X G C_5 and for every f 7^ 0. 

Lemma [3.21 is a modification of Lemma 7.2 of [23j, where the result is proved for p{x) = 
Pri\x\) and d = 3. 

Proof, (i) Let £x{x) be the fundamental solution of the operator —A + A^. In particular, if 

A I X I ~ 

d = 3, S\{x) = ■ We represent the entries of the matrix -D(A) as 




D^j{X)= / V,p{y)£x{y-z)Vjp{z)dydz. (3.11) 



The r.h.s. of (13. lip is defined and analytic for all A G C, then Dij{X) is defined and has an 
analytic continuation to C. Therefore, property (i) follows. 

(ii) By (13. lip . Dij{X) — > as |A| — )■ 00 with A G C_^. This implies property (ii). 

(iii) Note first that the matrix A{X) is positive definite for ImA = 0. Indeed, by condi- 
tion Al, we have 

V ■ A{X)v = iX' + u')\v\' - / >co'\v\'-v. Kov > 0, 



for any v ^ 0. Second, for ?/ G \ {0}, we find 



DiAiy + = - — / ^'^^ ^' = / ^^^^^ dr, 



where 

|2 



\e\=i 

\N\ 



Note that gij G C°°([0, +C)o)), and max (1 + r) \gij{r)\ < 00 for any > by condition A2. 

re[0,oo) 

Applying the Plemelj formula (see [15]) and condition A3, we obtain 

vlmDity + 0)v = -^y'i2n)-' J (v 9Y\p{\y\e)\^ dSe ^ (3.12) 

\9\ = 1 
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for any v e \ {0}. Hence, v ■ Im A{iy + 0)f = —v ■ lmD{iy + 0)^; 7^ 0. Therefore, there exists 
a (5 > such that v ■ lmA{iy + x)v 7^ for |x| < 6. Lemma [3.21 is proved. ■ 

Denote by A/'(t) the inverse Laplace transformation of A/'(A), 



e^*A/'(A) dX for t > 0. 



(3.13) 



~ joo— <5 



Lemma 3.3 Let m = 0, and let conditions A1-A3 hold. Then, for j = 0, 1, 

\^^'\t)\ < Ce-^\ t>l. 



(3.14) 



Proof. For j = 0, the bound on Af{t) follows by Lemma [3.21 To prove the bound for Af{t), 
we consider XJ\f{X) and prove the bound 



V ■ (W(A))'t; 



< for A G C_5. 

1 + Ar 



(3.15) 



Therefore, 



t^f{t) 



C 



e^*(AAr(A))'rfA 



acA=-(5 



-5t 



and bound (13.141) for Nit) follows. By Lemma [3^ (ii), to prove bound (I3.15P it suffices to 
show that < C{1 + |A|)-^ By formula fl3ll]) . we have 

D[A\)-^Q as |A| ^ cx) with A G C_5. 



Therefore, fl3.10p and Lemma [3.21 yield 

C 



1+ |A| 



+ |A| 



as lAI — )• 00. 



This implies (I3.15p . Bound (I3.14p with j > 2 is proved in a similar way. ■ 

Corollary 3.4 The solution of the Cauchy problem ( [J.^[ )-( rO]) has a form qt = Af{t)po + 
J\f{t)qo- Hence, Theorem \3.1\ follows from Lemma \3. 3[ if m = 0, and from Corollary \9.8\ if 
m 7^ 0. 



Corollary 3.5 Denote by V{t) a solving operator to the Cauchy problem ( 3.4)- (37E\) . Then 
the Duhamel representation for the solution of problem ^3. ( 15*. 5\) gives 



Qt 



V{t) 



Po 



+ I V{s) 



and F{t) = (Vpo; W/j0o)- Hence, Theorem \3. 1\ implies the bound for the solutions of l\3.^) : 



\qt\ + \qt\ < Ci£„(t)(|go| + bo|) + C2 j s^is)\F{t - s)\ ds. 





(3.17) 
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3.2 Time decay of Yt when (l)o{x) = for > Rq 

For the solution Yt to {\2.3\i . the following bound holds. 

Lemma 3.6 Let conditions A1-A3 hold and let Yq E S be such that 

fo{x) = 7ro(x) = for \x\ > Rq, (3.18) 
with some Ro > 0. Then there exists a constant C = C{R, Rq) > such that 

\\Yt\\£,R<CeUt)\\Yo\\£,R,, t>0, (3.19) 
for every R > 0, where Emit) is defined in /1 3. 6\) . 

Proof Step (i): At first, we prove bound (13.191) for '■= {<lt,Pt)- Let m = 0. Then, by 
condition (I3.18P and the strong Huygens principle, 

Wtcpo = for xe Br and t> R + Rq. (3.20) 

Since F{t) = (Vpo,W^i0o) = -(Po, VW^t^o), bound iKTTh and condition A2 imply 



\^t\ < C,e~''\^o\ + C{p)Je~'^\\V{Wt-s<l>onL^iB.)ds 





t 



< C,e-''\^o\ + C2 J e-''\\MRods<Ce-''\\Yo\\e,R, 

t~Ro~Rp 



(3.21) 



for t > Rq + Rp. Let m 7^ 0. Instead of (13.201) . we apply the well-known bound (rf = 3): 

\mMR<C{l + t)-'/^MRo, t>0. (3.22) 
Hence, \F{t)\ < C{1 + t)-3/2||0p||^^^ and bound fl3J9|) for follows from 

Step (ii): Now we prove bound (13.191) for (pf. Let m = 0. Then Eqn. (13.11) . condition (I3.18p . 
bounds (13^01) and (Km yield 

t 

Ut\\R<C J e-''\\Yo\\e,R,ds<Ce-''\\Yo\\e,R, for t > R + max{Rp, Ro} . 

t-Rp-R 

Let m 7^ 0. Therefore, Eqn. (13. ip . bound (I3.22p . estimate (13.190 for qt yield the inequality 
WMr < C{l + t)-^/^\\Yo\\£,Ro. This proves LemmaEH ■ 

4 Compactness of measures fit 

Proposition 12.111 can be deduced from bound (14.11) below by the Prokhorov Theorem [29| 
Lemma II. 3.1] using the method of [291 Theorem XII. 5. 2], since the embedding S = £^ G 
is compact for every £ > 0. 

Lemma 4.1 Let conditions A1-A3 and S0-S2 hold. Then 

supE||StFo||f < C'(^) < 00, Vi?>0. (4.1) 
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Proof. Let us write = e^°*, where Cq is defined in f l8.3p . At first, note that 



supE\\S^Yo\\iji<CiR), Vi?>0. 

t>o 



(4.2) 



Indeed, by the notation f l2.5p . II^^Yolll,/? = llW^t^ollfi + kt'P + where g° is a solution to 
the Cauchy problem 

g° + a;2g0 ^ g, t e M, (g?, g?)|t=o = (go,Po). 
Hence, + < C(|go| + \Po\)- By [F, Proposition 3.2] and [TOl Proposition 3.1], we have 

supE\\Wt(j)o\\l<C{R), Vi?>0. (4.3) 

t>o 

It implies (14. 2p . Further, we represent the solution to (12. 3p as 

t 

StYo = + j St^sBS%ds, (4.4) 



where the operator B is defined in (18. 3p . Hence, (I3.19P and (14. 2 p yield 

E\\StYo\\l^ < E\\SfYo\\ln + Ej\\St^,BS%\\lj,ds 



< C{R) + }sl{t-s)E\\S%\\l^^ds<C,{R)<^. U 



5 Asymptotic behavior of Yt in mean 

Proposition 5.1 Let d = 3, and let conditions A1-A3 and S0-S3 be satisfied. Then 

(i) the following hounds hold, 

E\qt-{Wt(t)o,a)\^ < Csmit), (5.1) 
Ebt-(Wi0o,/3)|' < CeUt), t>l, (5.2) 

where Emit) = e"^*^* if m = 0, and Emit) = (1 + t)~^ otherwise, the functions a and (3 are 
defined in Ii2.18\) and Ii2.19\) . 

(ii) Let f e [C^{R^)]'^ with supp / C Br. Then, for t > 1, 



E 



bt,f)-{Wt<Po,f*) <Cem{t) 



(5.3) 



where the function f^ is defined in Ii2.20\) . 

Proof, (i) At first, representation (I3.16p . Theorem 13.11 yield 



E 



qt- J A/'(s) {Wt-s(l>o, Vpo) ds < Celit) 





(5.4) 
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with Emit) from (13. 6p . Further, 



E 



+00 



+00 



Mr{Sl)dSi y"A/;r(s2)E(^(Wt_,,0o, VrPo)(W^t-S20O, VrPo)) 



For any t, Si,S2 G M, we have 

E({Wt.s,(pO,'^rPo){Wt-s,^0,'^rPo)) < C SUp E| (1^,00, V,Po) |' 

< CiSUpEWWrMR < C2 < OO 

by bound (I4.3p . Hence, applying bounds (13 .14^ and ( 19.12P for J\f{t), we obtain 

2 / I ^ \ 2 



E 



+00 



Af{s){Wt-s(l)o,'^Po)ds 



< / £m(s) ds = Ci£„(t). 



Therefore, (15. ip follows from (15. 4p . (15. 5p and (I2.18P because 

{Wt.^c^o, Vpo) = {Wtc^o, W^isVpo) • 
The bound (15. 2p can be proved in a similar way. 

(ii) Let / e [C^{R^)f with supp/ C Br. By Eqn. CT . 

t 

{cl^t, f) = {Wt^o, f)+ [ qt-s ■ (W^.Vp°, /> ds. 



for s> Rp + R if m = 0, 
+ s)-3/2) if m 7^ 0. 



E 



y (qt^s - (W^t-.0o,a)) ■ (W^sVp°,/> ds 



< 



(5.5) 



(5.6) 



(5.7) 



Note that 

{WsVp\ f) = 

Applying the Parseval inequality and bounds (15. ip and (15.70 . we obtain 

y\eUt-s)y/'\{W,Vp',f)\ds 
Jo 



Indeed, for m 7^ 0, by direct calculation, 

[\eUt - s)y/'\ (iy.Vp°, f)\ds<C [\l + t - s)-'/\l + sy/' ds <C,{1 + t)-^/\ 
Jo Jo 

00 2 

Write I{t) := E / {Wt^s(po, o) ■ (iy,Vp°, f)ds . If m = 0, then I{t) = for t > Rp + R hj 

t 

fl5TD . If m 7^ 0, then 



(5.^ 



\Ht)\<CEUt)- 



(5.9) 
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This follows from (15. 7p and from the following estimate: 



3 3 

^{WM,a)\^ = 5^ / Mir{s) {Wr-s(i>0,VrPo) ds 



i=l r=l 



< C < oo, for r G 



by dH, (EH and (ITOl) . Relation (EE]) and bounds (ESD and (ES]) imply ( lOll . 
Corollary 5.2 Let Z = (/, m, t;) G P = [Co°°(M3)]2 x x M^, y/^g^ 

{Yt,Z) = {WtcP,MZ))+r{t), 



where Vl{Z) is defined in \2.11/^ , (F^, Z) = (0^, f) + qfU+pfV, Yt = is a solution 

to the Cauchy problem / fO|) . and E (|r(t)p) < CSmit). 



6 Convergence of characteristic functionals 

Proof of Proposition 12.121 By the triangle inequality, 



< 



^^i{Wt<i>o,n{z)) 



e 2- 



.(Z,Z) 



The first term in the r.h.s. of (16. ip is estimated by 



E 



{Yt,Z) - {Wt(l)o,niZ)) < E\r{t)\ < (E|r(t)|2j < CeH^it) ^ as t 



— )■ oo, 



■ (6.1^ 



(6.2) 



by Corollary 15.21 It remains to prove the convergence Kexp{i{Wt4>o,^{Z))} = fifiVL^Z)) to a 
limit as t — > oo. 

In [HI [in], we have proved the convergence of fif{f) to a limit for / e = [C^(M'^)]^. 
However, VL{Z) ^ Vq in general. Now we introduce a space H'^ such that fi(Z) G H'^ and the 
characteristic functionals fif{f), t G M, are equicontinuous in H^. As before, ^(A;) stands for 

(1^1 



2N1/2 



Definition 6.1 H!^ is the space of the pairs f = (/", /^) of real-valued functions /° and f^ 
such that /°, e L2(]R3)^ ^ ifl(]R3)^ ^^^^/^ ^;^g y^^^^g ^^^^ 



2 ._ 



\\fr+\\-^-\k)f%k)r+\\fr+\\vfr- 



If m ^ 0, H'^ = L2(R3) © H\M.^). Note that fi(Z) G H'^ if Z e V. This follows from 
formulas (I2.17p - (l2.20p and from condition A2. 

Remark 6.2 By conditions SO, S2 and S3, the correlation functions Qq (x,?/) of the measure 



Hq satisfy the following bound: for a, /3 G , \a\ < 1 



< 1 - J, hj = 0,1, 



\D:i^Q'^ix,y)\ < Ceoif'/'ilx -y\), x,y e R^ 
according to [TSl Lemma 17.2.3]. Therefore, by (12. 9p . 

oo 

D:iQ^^{x,y)\^dy<Cel [ ^^/\\x - y\) dy < C^e^ [ rV^' {r)dr < oc, p>l 



(6.3) 



Hence, by the Shur lemma, the quadratic form {Qo{x, y), f{x)^f{y)) is continuous in [L^(R'^)]^. 
Moreover, the quadratic form Q^{f,f) is also continuous in [L^(R'^)]^. 
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Lemma 6.3 (i) The quadratic form Qf{f,f) := j \{(f)o, f)\'^iJ,f{d(f)o), t ^M., are equicontinu- 
ous on H'^. (a) The characteristic functionals /tf (/), t G M, are equicontinuous on H'^. 

Proof, (i) It suffices to prove the uniform bound 

sup|Qf(/,/)|<C||/||L feH'^. (6.4) 

Note that Qf{fJ) = {Qo{x,y),Wlf{x)®Wlf{y)). Then, by RemarkE^l 



sup|Qf(/,/)|<Csup||iy;/||i. <C||/|| 



I 



since in the Fourier transform, we have 
(ii) By the Cauchy-Schwartz inequahty and (16. 4p . we obtain 



cosLo{k)t -uj{k)smLo{k)t \ f f°{k) 
oo-\k) sin 00 {k)t cosuj{k)t J \ f\k) 



for / = (r,r). (6.5) 



Af(/i)-/if(/2) 



/ii {d(f)o) 



< / 1(00, /l - /2)|/^f(#o) < W / 1(00, /l ~ /2)IX(rf0o) 



\/2f(/l-/2,/l-/2) < C\\h - f. 



2 m- 



We return to the proof of Proposition 12.121 Since Q{Z) G H'^, Proposition 2.3 in [8] (or 
Proposition 3.3 in ^Q\) and Lemma [6.31 (ii) yield 

fif{n{Z))^exp{~QZ{n{Z),n{Z))} as t^oo, 

where is defined in (12.161) . This completes the proof of Proposition 12. 121 and Theorem 12. 1U[ 



6.1 Convergence of correlation functions 

Proposition 6.4 Let all assumptions of Theorem \2.10\ be fulfilled. Then, for Zi, Z2 G V, 

E {{Yt, Z,){Yt, Z2)) ^ Qoo(^i, Z2), t^ 00, (6.6) 
where Qoo{Z, Z) is defind in IH 



Proof. It suffices to prove the convergence of EKy^, Z)\^ to a limit as t — )■ 00. It follows from 
Corollary [El that for Z eV, 

n{Yt, Z)\^ = n{Wt(t>o, ^{Z))\^ + 0(1) = Qf (l^(Z), n{Z)) + 0(1), t ^ 00, 



where Vt{Z) is defined in (I2.17P and Vt{Z) G H'^^. Therefore, by Lemma 5.1 in [8] (or Propo- 
sition 6.1 in [lO]) and by Lemma [631 (i), we have lim Qf (1](Z), fi(Z)) = Q^(fi(Z), fi(Z)). 

Formula ^TM imply M . 
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7 Ergodicity and mixing for limit measures 

Denote by Eqo (E^) the integral w.r.t. fioo (z^,^, respectively). In [5J, we proved that Wt is 
mixing w.r.t. /x^, i.e., for any /, (7 G L2{'H, fi^), the following convergence holds, 

EZifiWt<l))gm^EgfEgg as t ^ 00. (7.1) 

Recall that the limit measure /loo is invariant by Theorem 12 . 1 01 (iv) . Now we prove the mixing 
property for the flow St- 

Theorem 7.1 The flow St is mixing w.r.t. fioo, i-^-, for any functions f,g & L2{S,fioo), 
have 

Eoo {f{StY)giY)) ^ Eoo/ E^g as t ^ 00. 
In particular, the flow St is ergodic w.r.t. /loo, i.e., for any / G L2{£,fioo), 

T 

\im ^ [ f{StY)dt = E^f{Y) (mod/ioo). 

T— >CXD I J 



To prove Theorem 17. H we introduce new notations. 

Represent Y e £ &s Y = (F°,r^) with Y^ = (v?(-),9) e Hl^(R'^) x R^, Y^ = (7r(-),j9) G 
Ll^iR"^) X R'^, andZ eV as Z = {Z°, Z^) with Z^ = (/O(-), u°), Z^ = {f\-), u^) G C^{R'^) x 
M*^. For t G M, introduce a "formal adjoint" operator S'^ on the space "D, by the rule 

{StY,Z) = {Y,S',Z), YeS, ZeV. (7.2) 



Lemma 7.2 For Z eV, 

S',Z={fti-),UtJti-),Ut), (7.3) 
where {ft{x),Ut) is the solution of system h2.1\) with the initial data (see 112.^) ) {fo,qo,'^o,Po) = 

(r,n\/°,n°). 

Proof Differentiating in t with Y,Z eV, we obtain (S'tF, Z) = (F, ^^Z). The group St 
has the generator 

The generator of 5*^ is the conjugate operator C = ^ ^ ) Hence, (17. 3p holds with 



Ut J \ Ut 

Since the limit measure /Xqo is Gaussian with zero mean, the proof of Theorem 17.11 reduces 
to a proof of the following convergence. 

Lemma 7.3 For any ^1, ^'2 G T), 

Eoo((^tF,Zi)(y,Z2)) ^0, t^oo. (7.5) 
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Proof First we note that, by relation ( I2.22p . 

((F, z,) (F, z^)) = eZ ((0, f^(^i)) (0, niz,)) 

where Q{Z) is defined in (12.171) . Secondly, for fixed t, we have S'f-Z G P. Further 
E^(^{StY,Z^){Y,Z2)) = E^, ((F, 5,%) (F, Z2)) = ((0, (0, r](Z2)) 

= E^f(0, ins[-win)z^){<p,niZ2))) +Eg({<p,w;niz,)){<p,niz,)) 



= h + h. (7.6) 
Note that (0,fi(^)) G L2(H,/i^) for all Z G P. Indeed, by fl2:22D . 

E^|(0,n(Z))|2 = {Q'Z,n{Z)®n{Z)) = {Q^,Z®Z) < 00. 
Therefore, the convergence fl7.ip implies 

h ^ E^((0,iy;r](Zi))(0,i](Z2))^ 

= Ef,f(iy,0,fi(Zi))(0,fi(Z2))) ^E^(0,n(Zi))E^(0,n(Z2)), t^oo. 



On the other hand, 

E^(0,fi(Z,)) =E„o(F,Z,) = 0, for Gl), 
because /Xoo has zero mean. Therefore, 

/2 -> 0, t ^ 00. (7.7) 

Now we prove that E^|(0, Vt{S[Z) - W/fi(Z))|2 = for all t > 0. This yields 

h ^ E^((0, {ns[ - w;n)z^){<p,n{Z2))) = o. (7.8) 



Indeed, by Corollary 15.21 

E\{Sr+tY,Z) - {Wr+t<P,^{Z))\^ ^0, r ^ oo. 
On the other hand, since {Sr+tY, Z) = {SrY, S[Z), we have, for all t > 0, 

E| {SrY, S[Z) - {Wr(j), ^{SlZ))]"^ -> 0, r ^ oo. 
Therefore, by the triangle inequality, 

A := E\{Wr(j),n{SlZ)) - (l^,+t0,fi(Z))|2 ^0, T^oo. 
Since, (VT^+t 0, ^](Z)) = (W^, 0, iy;fi(Z)), we obtain 

A = E\{Wr(j), ^{S'tZ) - Wln{Z))\'^ ^0, T^oo. 
By Theorem 12.81 (v) and Lemma 16.31 

E^l {(f), n{StZ) - w;n{z))\'^ = lim E\{Wr(j),n{StZ) - w;n{z))\'^ = o for aii t > o. 

T—^OO 

Finally, fl7.6p -( rrH]) imply the convergence (17. 5p and then. Theorem 17.11 
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8 Appendix A: Existence of solutions 

In this section, we prove Proposition 12.21 by using the methods of [231 Lemma 6.3]. 

Represent the solution Yt as the pair of the functions (Y^,Y^^), where Y^^ = {ft, It), Yt - 

{T^t,Pt)- 

Let E be the Hilbert space of pairs Y = where F° = {(p,q) G H^M:^) © 

llV^f + m2||^|p + ||vr|p + |g|2 + |pp 



Now we prove the axillary lemma 



with the finite norm : = 



Lemma 8.1 Let conditions Al' and A2 hold. Then 

(i) for every Yq G E, the Cauchy problem 1^2. 3\) has a unique solution Yt G C(M, i?). 
(a) For every t G M, the operator 5*^:1^0'"^ is continuous on E. 
(Hi) The energy is conserved and finite, 

H(Yt) = H{Yo) for t G M. 



U) 



Proof. Step (i) li p = 0, then the existence and uniqueness of the solution Yt G C(R,E) 
to problem (12. 3p is proved by using the Fourier transform. Therefore, problem (12. 3p for Yt G 
C(M, E) is equivalent to 



Yt = 6^°% + / e^'^'^'^BYs ds, 



where (cf fTTiD ) 

^0 = 



B{Y'^,Y^' 



I 

-Ao 

= (0, RY 



A 



q 

RY^ : 



—A + m'^)(p 
q ■ Vp, {ifi, Vp) 



for F° = {ip,q). 



1.3) 



Note that ||e^"*yo||£; < Cll^olU; and the second term in (18. 2p is estimated by 

t 

sup II / e^''>^'-''^BYsds\\E < CT sup ||F,| 

\t\<T J \s\<T 



s\\E- 



This bound and the contraction mapping principle imply the existence and uniqueness of the 
local solution Yt G C{\—e,e\,E) for some e > 0. 

Step (a) To prove flHTTl) . we first assume that 0o = {vo, ^o) e C^(M^) xC^(]R^) and 0o(a;) = 
for |x| > i?o- Then ^t{x) G C^i^l x M^) and 



for |x| > \t\ + max{_Ro, -Rp} 



by the integral representation (18.21) and condition A2. Therefore, for such initial data, relation 
( 18. ip can be proved by integrating by parts. Hence, the energy conservation ( 18. ip follows from 
the continuity of St and from the fact that CliM?) © © C'ol^^) ® is dense in E. 

Step (Hi) Note that 



^l|V^(x)f + ^m2||y,(x)f + (g-V^,p) 



2(27r)s 



+ m?Lp{k) + 



iq ■ kp{k) 



q ■ kp{k) 
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Hence, the Hamiltonian functional H{Y) is nonnegative. Indeed, 



^l^') = -{llvrf + (27r)-=^||VF + m2^(fc) + 



TFT 



+^(c^'|g|'-g-^mg) + ^br>0 (8.4) 

by condition Al'. Moreover, by flS.ip and fl8.4p . for |t| < e, we obtain 

m\l<CH{Yt) = CH{Y^). 
On the other hand, since \q ■ {V(p,p)\ < |(||Vv?P + I^zPHpIP), 

H{Y) < II V^f + ^m^ll^f + Iw^r + lico' + llpDigr + 1\P\'< C\\Y\\l, 
where C = max{l,m2/2, (o;^ + ||p|p)/2}. Hence, we obtain the a priori estimate 

\\Yt\\E < CilllolU for t e M. (8.5) 

Therefore, properties (i)-(iii) of Lemma [8. II for arbitrary t G M follow from bound (18.51) . ■ 

We return to the proof of Proposition 12.21 Let us choose R > Rp with Rp from condition 
A2. Then, by the integral representation (18.21) . the solution Yt for |x| < R depends only on the 
initial data Yq{x) with |x| < R + \t\. Thus, the continuity of St in S follows from the continuity 
in E. 

For every R> 0, define the local energy seminorms by 



1^ \\e{r) ■- 



jV<^(x)|2 + m^\if{x)\'' + |7r(x)|2) dx + \q\^ + \p\^, Y = {ip, vr, q,p). 

\x\<R 

By estimate (18.51) . we obtain the following local energy estimates: 

< C||Fo|||(fl+|t|) for R>Rp and teR. (8.6) 
Hence, bound (12. 6 p follows from (18. 6p . ■ 
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9 Appendix B: Time decay of qt and pt in the case m 7^ 

Theorem 9.1 Let m ^ and d = 3, let the conditions A1-A3 hold, and let Yq E E and 
Ii3.18\) hold. Then qt and pt are continuous, and the following bound holds: 

+ <C(l + t)-3/2, t>0. (9.1) 

This theorem is a modification of Lemma 17.1 in [19] . 

Denote by gx the fundamental solution of the operator —A + + A^, 

p-«;|?/l 

where = + A^, Re/t > for ReA > 0. Applying the Laplace transform, by ( ]3.4p -( !375|) . 
we obtain 



M(A) 
where 



and the matrix -D(A) is defined in (13. 8p . The entries of -D(A) are of the form 
The following result was proved in p!^, p. 351]. 

Lemma 9.2 (i) The operator —A + + A^ is invertible in L^(M^) for ReA > and its 
fundamental solution Ii9.^) decays exponentially as \y\ -> oo. 

(a) For any fixed y ^ 0, the Green function g\{y) admits an analytic continuation (w.r.t. 
the variable X) to the Riemann surface of the algebraic function ^/X^ + m? with the branching 
points ±im ifrriy^O. 

Lemma 19.21 and (19.31) imply that M(A) admits an analytic continuation from the domain 
Re A > to the Riemann surface with the branching points ±im with m ^ 0. Moreover, the 
matrix M~^(A) exists for large A. It follows from (19. 3p since Dij{X) — )■ as Re A ^ cxd by (19. 4p . 

Corollary 9.3 (i) The matrix M(A) is invertible for Re A > 0, and 

(i(A))-^^-'w(:)- ^^'>°- 

(a) The matrix M~^{X) admits meromorphic continuation from the domain ReA > to the 
Riemann surface with the branching points ±im with m ^ 0. 

Now we investigate the limit values of M~^(A) at the imaginary axis X = ix, x E M.. The 
limit matrix 

exists, and its entries are continuous functions of x G M, smooth for |x| < m, and for |x| > m. 
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Lemma 9.4 The limit matrix M{ix + 0) is invertible for x e M. 

Proof, (i) Let < m. Then the matrix (u;^ — x'^)I — D{ix + 0) is positive definite. Indeed, 
by condition Al, 

v.ii.^-xy-Di^x + 0))v = (^'-^M-i:^-)-'J ^"j^'^l^ff 

> {co'^ - m^)\vf - V ■ Kqv > (9.7) 

for every v e \ {0}. Therefore, M{ix + 0)^ 7^ for all ^ G \ {0}. 

(ii) Let |a;| > m. In this case. Then v • lmD(ix + 0)v 7^ for every e \ {0}. Indeed, 

-3 f kikj\f){k)\^ dk 



lmbij{ix + 0) = Im (27r)-^ j 
For £ > 0, consider the function 



kP' + rn^ — {x — iO)" 



d{ix + e) := /' \ . dk, \x\>m. 

' k-^ -\-m^ — [x — ley 



Denote h^{k) — A;^ + m^ — (x — is)^. For |a;| > m, 6o(^) = if \k\ = y/x'^ — vn?. Let us fix a small 
5 > and introduce a cutoff function C e C^{^^) such that C{k) > 0, C(^) = 1 if \bo{k)\ < S 
and = if |&o(^)| > 25. We split the integral d{ix + e) into two terms ds{ix + e) and 
fs{ix + e): 

d{tx + e)^ ! C{k)^^^^^ dk + /(I - C(fc)) ^-^^'f^f ' dk = ds{ix + e) + U{ix + e). 

J Oe{k) J he{k) 

Note that for fixed 5 > 0, 
f,ii:, + e):^J{l-C{k))^^^^^dk^ J {1 - ({k))^^^^^ dk e R as e ^ 0. 

\bo{k)\>S 

Therefore, Imd(ia; + 0) = lmds{ix + 0), where 

ds(ix + 0) = lim ds(ix + e), dgHx + e) :^ f CW ^'^/'ff dk. 

e^+O J he[k) 

Denote a{k) — y/k'^ + m?- Then h^{k) — a?{k) — {x — iey. Let x > 0. Since 

1 1 1 

+ 



b,{k) 2a{k){a{k) - X + ie) 2a{k){a{k) + x - ie)' 
then lmds{ix + 0) = lmd^{ix + 0), where, by definition. 

We rewrite di{ix + e) as 



a(fe)— a;=a 
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Hence, lmd^{ix + 0) = —TTg{0) by the Plemelj formula. Finally, for x > m > 0, 

Therefore, applying condition A3 we obtain, for |x| > m, 

T;-ImD(zx + 0)t; = -sign(a:)7r(27r)-3 I rfS ^ (9.8) 

J 2|A;| 

\k\=\/x'^—m^ 

by condition A3. Thus, formula (19. 6p and estimate (19. 8p imply Lemma 19.41 Indeed, assume 
that there exists a nonzero vector v = (a,/3) G M'^ x such that M{ix + 0)v = 0. Then, by 
(19. 6p . we obtain 

ixa — f3 = \ P ~ 

[u'^I - b{ix + Q)^a + ix(5 = Q ^ j {{u'^ - x"^)! - D{ix + Q)^ a = 

Hence, ([92D and ([ISD imply that a = 0, and then, /3 = 0. Therefore, KerM(zx + 0) = 0. 
Lemma 19.41 is proved. ■ 

Remark 9.5 (i) In the case when m ^ condition A2 can be weakened. It suffices to assume 

the real-valued function p to be in the Sobolev space H"^ and of compact support. 

(a) We use condition A3 only in estimates 1^3.12) and 1^9. 8\) . Hence, instead of condition A3, 

it suffices to assume that, for any v eM.^\ {0} and for any r > 0, / (f ■ 9)'^\p{r 9)\'^ dSg ^ 0. 

|0|=i 

Corollary 9.6 The matrix M~^[ix + 0) is smooth w.r.t. a; G M outside the points x = ±im, 
with m 7^ 0. 

Lemma 9.7 (i) The matrix M~^{ix + 0) admits the following Puiseux expansion in a neigh- 
borhood of ±im (m 7^ 0): there exists an e± > such that 

oo 

M-\ix + 0) = '^c^{xTm)''/'^, \xTm\<^±, a; G M. (9.9) 

(a) There exist a matrix Rq and a matrix-valued function Ri{x) such that 

M-\ix + 0) = -Ro + Ri{x), |x|>m + l, x G M, (9.10) 

X 

where \d^Ri{x)\ < Ck/\x\'^ for \x\ > m + 1, x e R, k = 0,1,. . .. 

Proof, (i) Formula ([91]) follows from and flQ]) . 

(ii) Let / G L\R^) with supp / C Bn- Then (see formula (16.7) in [l9]) 

\\d'^[-A + m' + {zx + 0)r'fL.,s.^ < ^^ll/IUw |x| > m + 1, 



x\ 



for every R> 0. Therefore (see (13.81) ) we obtain \d^D{ix -\- 0)\ < Ck/\x\ for |x| > m + 1. Then 
formula 1^ implies i^M>- ■ 
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Proof of Theorem 19.11 Applying ( I9.5p . we obtain 

+00 

*J._L|,-M-(» + o)(* (9.11; 

—00 

We split the Fourier integral f l9.1ip into two terms, by using the partition of unity: Co (2;) + 
Ci(x) = 1, X e M: 

+00 

1)=^ I (Co(^) + Ci M-\^x + 0) ) dx = loit) + hit), 

—00 

where the functions Cfc(^) ^ C°°(]R) and suppCi C {s G M : m/2 < |x| < m + 1}. Then (i) 
the function Io{t) G C[0, +00) decays faster than any power of t due to Lemma [9.71 and bound 
fl9Tl) . and (ii) the function Ii{t) G C°^(M) decay like (1 + \t\)-^/^ by dH]). Theorem EH is 
proved. ■ 

Corollary 9.8 Let m 7^ and conditions A1-A3 hold. Then 

\Af^'\t)\<C{l + \t\)-'/', J = 0,1, (9.12) 

where Af{t) is the inverse Laplace transform of J\f{X). This bound can be proved in a similar 
way as Theorem \9.1[ 
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10 Appendix C: Non translation invariant initial mea- 
sures 

The proof of Theorem 12.101 is divided into two main steps: first, we prove the asymptotic 
behavior of solutions Yt in mean: {Yt, Z) ~ {Wt(f),^l{Z)) as t — )■ oo (see Corollary 15. 2p . This 
asymptotics allows us to reduce the convergence analysis for the coupled system to the same 
problem for the free wave (or Klein-Gordon) equation. Second, we prove the weak convergence 
of the measures /if = W^fiQ as t — )■ oo (see Theorem 12. Sp . 

However, the weak convergence of /xf holds under weaker conditions than S2 and S3. Now 
we formulate these conditions on Hq (see [8] if m = and [ID] if m 7^ 0). 



10.1 Conditions on /Iq 

In the case m 7^ 0, we assume that /x^ has zero mean, satisfies a mixing condition S3 and has 
a finite mean energy density (see (12. 7p ). i.e., 

J (|¥Po(x)r + |Vv^o(x)r + |vro(x)|2)/io^(d0o) = g[J°(x, x) + [V,V,Q°°(x, y)] [^^ + Qji(x, x) 

< eo < 00. (10.1) 

However, condition S2 (see Section 12. 3p can be weakened as follows. 
S2' The correlation functions of the measure fi^ have the form 

Qi/(x, y) = <il{x - y)C-(^i)C-(yi) + - 2/)C+(^i)C+(2/i), hi = 0, 1. (10.2) 

Here q±{x — y) are the correlation matrices of some translation-invariant measures fi^ with 
zero mean value in x = (xi, . . . , x^), y = {yi, ■ ■ ■ ,yd) ^ a > 0. The measure /Iq is not 
translation- invariant if g!f 7^ 5'+. The functions (± G C°°(R) such that 

a(.) = (i' f (10.3) 

^ [ 0, for ± s < -a. ^ ^ 

If m = 0, instead of S2 we impose condition S2'. 
S2' The correlation functions of fi^ have the form 

with q± as in (110.20 . However, if m = 0, instead of (110. ip we impose a stronger condition 
(110. 5p . Namely, the following derivatives are continuous and the bounds hold. 



\D::SQi\x,y)\< 



Cu,i\x-y\) if /€ = 0,l,...,d-2 
Cud^idx — y\) if K = d—l,d,d+l 



K = i + j + \a\ + \(3\, (10.5) 



with |a| < (rf- 3)/2 + i, \f3\ < (d- 3)/2 + j, i,j = 0,1. Here z/^ G C[0,oo) {k = 0, . . . ,d - 1) 
denote some continuous nonnegative nonincreasing functions in [0, 00) with the finite integrals. 



1 + r)'^ ^Vi^{r)dr < 00. Moreover, for (i > 5, / (1 + rY ^^'^i'^{r)dr < 00 with k = 0,2 
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Now we formulate a new mixing condition S3' for the measure /z^ (if m = 0). 

S3' Let a = (ai,a2,«3) with integers ctj > 0. Denote by <7ia{A) the ci-algebra of the subsets 
in "H generated by all linear functionals 

0o = (0°,0j)H^(I)>^,/), with |a|<l-^, z = 0,l, 

where / G C^(M^) with supp / C A. For k = 0, 1, let o"k(^) be the cr- algebra generated by 
aia{A) with i + \a\ > K, i.e., af^{A) = V aia{A). 

i+\a\>K 

We define the (Ibragimov) mixing coefficient of a probability measure fi^ on "H as (cf (12.81) ) 

\fi^{AnB)-f,^{A)f,^{B)\ 

^^i,K2{r) = sup sup ^— , Ki,K2=0,l. 

f^l^iB) > 

We assume that the measure Hq satisfies the strong uniform Ibragimov mixing condition, i.e., 
for any Ki, K2 = 0, 1, V^ki,k2(^) — 0, r — )• oo. Moreover, 

'fKi,K2ir) < Cvl{r), where k = ki + Kg, Ki, ^2 = 0, 1. 



Remark 10.1 (i) In [HI [TO], we have constructed the generic examples of the initial measures 
IJjQ satisfying all assumptions imposed. 

(n) Condition S3 and the bound (110.11) imply the bound (16.31) . 

(iii) Condition (110.51) imphes (110. ip . Condition S3' implies estimates (110. 5p with i + |a| < 1, 

J + l/31 <1. 

(iv) The conditions S2 and S3 admit various modifications. We choose the variant which 
allows an application to the case of the Gibbs measures (see Section flLSp . This is related to 
the slow long-range decay of the correlation function Q^{x^y) ^ \x — yp"'^, — 1/| — ^ oo. 

The mixing condition S3' is weaker than condition S3. On the other hand, the estimates 
(110. 5p with K > 2 are not required for translation-invariant initial measures or in the case 
when m 7^ 0. 



10.2 Convergence to equilibrium 

Theorem 10.2 (see /3, \W^) Let all conditions imposed on fi^ hold. Then the assertions 
of Theorem \2.8\ remains true with the limit correlation matrix Q^{x,y) = q^{x — y) of the 
following form. In the Fourier transform, q^{k) = q^{k) + q^{k), where (cf l[2.1Ji\) ) 

q^ik) = i(q+(fc) + C'(A:)q+(A:)C7^(A:)), 

q^{k) = tsgn{k,)l(d{k)ci-{k)-ci-{k)C^{k)y 

with q"*" = |(g+ + g_), q^ := |(g+ — g_), and C{k) from h2.15\) . 

Corollary 10.3 Let conditions SO, SI and all assumptions of Theorem \10.2\ be fulfilled. Then 
the assertions of Theorem \2.1(A hold with introduced in Theorem \10.2[ 

In Section n 1.31 we give an application of Theorem 110.21 to the case when fi^ (see S2') are 
Gibbs measures with different temperatures T+ 7^ T_-. 
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11 Appendix D: Gibbs measures 
11.1 Phase space 

We define the weighted Sobolev spaces with any s, a G M. 

Definition 11.1 (i) if^(]R'^) is the Hilbert space of the distributions G S"(M'^) with finite 
norm 

= II (x)"AVII < oo, AV = [{kYip{k)] , s, a G M. (11.1) 

(a) l-il^ = H^^iW^) © H^^iW^) is the Hilbert space of pairs (p = {ip{x),n{x)) with real valued 
functions ip{x) and 7t{x), with finite norm 



a = llv^lls+i.Q + \H\s,a, s, a G M. (11.2) 

(Hi) = © R"' © M"' is the Hilbert space of vectors Y = {(f){x), q,p), with the norms 

\\y\\s,a = ||0||s,a + |g| + bl, S,aeR. 

Note that C 7/^ (and also £^ C if s > s and a > a, and this embedding is compact. 
Moreover, for any a, "H^ C Ti, £^ C S (see Definition 12. ip . 

Remark 11.2 Let s, o; G M, and conditions Al' and A2' hold. Then (i) for every Yq G 
the Cauchy problem h2.3\) has a unique solution Yf G C(M.,£^). 

(a) For every t G M, the operator St : Yq Yt is continuous on S'^. Moreover, there exist 
positive constants Ci,C2 > such that ||5'(FolUa < Ci(^)'"^ ll^olUa- 



This bound can be proved as in [2T], where the nonhnear "field-particle" system was studied. 
11.2 Gibbs measures for heat reservoir 

Write = (v2,vr). We introduce the (normalized) Gibbs measures on the space H^. For- 
mally, 



Now we adjust the definition of the Gibbs measures g^. Write = (0°,0^) = (v',7r), and 
denote by Q''^{x,y), i,j = 0, 1, the correlation functions of g^, 

Q''{x,y) = j 4>\xW{y)fif{d4>) = q'^{x-y), x,y eR''. 

We will define the Gibbs measures g^ as the Gaussian measures with the correlation functions 

q''\x-y)=TSrn{x-y), q'\x - y) = T6{x - y), q^\x - y) = q'\x - y) = 0, (11.3) 

where T = 1/(3, £m{x), x G M'^, is the fundamental solution of the operator —A + m^. The 
correlation functions g** do not satisfy condition (110. ip because of singularity at x = y. The 
singularity means that the measures g^ are not concentrated in the space "H. 
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/ G C^{W). (11.4) 



Definition 11.3 For /3 > 0, define the Gibbs measures g^{d(f)) as the Borel probability mea- 
sures g^{d(p) = g^pidif) x glidn) zn = H^+\R'^) (g) H^{R'^), s,a< -d/2, where gl{d^) and 
gjj{d7r) are Gaussian Borel probability measures in spaces H^~^^{M.'^) and H^{R'^), respectively, 
with characteristic Junctionals 

glU) = |exp{z(^,/)}^7°(rfy,) = exp|-^((-A + m2)-V,/)} 
glif) = y"exp{z(7r,/)}^?J(rf7r)=exp|-^(/,/)| 

By the Minlos theorem, the Borel probabihty measures g^ and gjj exist in the spaces 
if^+^(M'^) and H^{M.'^), respectively, because formally 

j M\l+i,a9p{d^) <oo, J WrrWl^ glidn) < oo, s,a<-d/2. (11.5) 

We verify ffTT3D . Definition ffTTTTD implies, for ip G H'^{R'^), 

||<^||2^ = (27r)-2'^ / {x)^"( [ e-'''^^-'''\kY{ky^{k)^{k')dkdk')dx. (11.6) 



Let g{dif) be a Gaussian translation invariant measure in H^{R'^) with a correlation function 
Q{x,y) = q{x — y). Let us introduce the following correlation function 



C{k,k')^ J m0{k')g{dip) 

in the sense of distributions. Since (p{x) is real-valued, we have 

C{k,k') = F,^kF,.^.u'Q{x,x') = {2TxY5{k - k')q{k). 
Then, integrating f lll.6p with respect to the measure g{d(p), we obtain the formula 

M'g{d^) = {27r)-''[ {xf-dx I {kY^q{k)dk. 



Substituting q{k) = T (see flll.Sp ) we obtain the second bound in flll.Sp . To obtain the first 
bound in flll.Sp we replace s into s + 1 and put q{k) = TS^nik) = T^kl"^ + m^)~^. 

Below for spaces £^ and 7^^, we put s,a < —d/2. 
Definition 111.31 implies the following lemma. 

Lemma 11.4 The Gibbs measures g^ are invariant w.r.t. Wt- Moreover, the flow Wt is 
mixing w.r.t. g^ . 

This lemma can be proved by the similar reasonings as in Sections 111.61 and [71 



28 



11.3 Application of Theorem 10.2 to Gibbs measures fi^ 

In in Ho], we studied the case when fi^ (see condition S2' and (110. ip ) are the Gibbs measures 
corresponding to different positive temperatures T_ "^T^. We define the Gibbs measures 
g^ as the Gaussian measures with the correlation functions (cf. flll.3p ) 

(^l{x-y)=T^Em{x-y). ql{x - y) =T^b{x - y), qf{x-y) = q'^{x-y) = Q, (11.7) 

where x,y E M."^. The Gibbs measures g^{d(f)) = g^^{d4>) (with I3± = l/T±) are the Borel 
probability measures on the space "H^, its definition is specified in Definition 111.31 

Let us introduce (0_, 0+) as a unit random function in the probability space ("H^ x "H^, g^ x 
g^). Then (f)± are Gaussian independent vectors in "H^. Define a Borel probability measure 
fiQ = g^ on "H^ as the distribution of the random function 

Mx) = C-(a;i)0_(x) + C+(a;i)0+(x), 

where functions (± are introduced in (110. Sp . Then correlation functions of gQ are of the form 
(110. 2p . where are the correlation functions defined in (lll.7p . Hence, the measure g^ has zero 
mean and satisfies condition (110. 2p . However, g^ does not satisfy (110. ip because of singularity 
a.t X = y. Therefore, Theorem 110.21 cannot be applied directly to /i^ = g^. The embedding 
"H^ C 'H'^ is continuous by the standard arguments of pseudo differential equations, [17]. The 
next lemma follows by Fourier transform and the finite speed of propagation for the wave or 
Klein-Gordon equation. 

Lemma 11.5 The operators Wt : (po ^ (j^t allow a continuous extension W ^ W . 

Let 00 be the random function with the distribution g^ . Hence 0o £ a.s. Denote by gf 
the distribution of PVt0o- 

Theorem 11.6 Let s < —d + 1/2. Then there exists a Gaussian Borel probability measure g^ 
on the space l-L^ such that 

g^^gg, t^oo. (11.8) 

Proof Let us fix an s < —d+1/2 and introduce the random function 0q := A'^^o- By definition, 
(see (dllD) 

A>o = As(-) * 00, where A,(a;) := F^^^{{ky). 
Since 0o ^ T^^, <Po ^ T^a ^-S- Let us denote by g^ the distribution of Wt{A'^(po), t gM., 

gf^'{uj):=gl^^'{Wi'uj)=gf{A-'uj), ioiueBiUl), t G M. 

Then gf'' = gfh'^ and gf = gf^'k^ since K^Wtcj),) = Wt{h''(t>o). 

Let Ql{x,y) stand for the (matrix) correlation function of the measure g^ 

Ql{x,y) = j (0^(x)®0g(y))(7f'^(rf0^) = j (0o(x)®0o(l/))^7f(^^0o)*(A.(x)A,(y)). 

Now we study properties of g^'^- The measure (^jf obviously has zero mean. The correlation 
matrix of g^'^ is of the form 

Q^(x, y) = Qo{x, y) * (A,(a;)A,(y)). (11.9) 
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The correlation function Ql{x,y) satisfies fll0.4l) with a suitable modification: fll0.4p holds up 
to (5(1 + |x| + \y\)~^ with any 5, > and with a = a{6, N). Namely, 



Qoix,y) 



Q-ix - y) + R-{x,y), xi,yi < -a 
ql{x - y) + R+{x,y), Xi,yi>a, 



where q^{x — y) := q±{x — y) * {As{x)As{y)) are correlation functions of the random function 
A^4'±, and \R±{x,y)\ < C(l + |x| + |?/|)~^ with any A^. This follows from the convolution 
representation flll.Qp . since As{x) is a function of class //[^^(R'^) with a rapid long-range decay. 

Now we verify condition (110. 5p . If m = 0, q^{x) = C±\x\'^~'^ (since d > 3). Therefore, 
condition (110. 5p also holds for g^'^ with the functions //^(r) = C(l -f t-^-'^+^-k for a sufficiently 
large C = C{s,T±) > 0. It follows immediately for s < with sufficiently large \s\ from the 
same convolution representation (lll.9p . For s < —d + 1/2 it follows by the pseudodifferential 
operators techniques. 

Since gf'^ is Gaussian, we do not need the mixing condition S3 or S3'. Applying the 
technique from |8l [10], we first prove that the correlation functions Qf{x,y) of the measures 
gf'"^ converge to a limit as t — )■ cxd. Further, we derive the following bound: 

for any R> 0, sup / ||0o||ii5'f ''*(c^0o) < C < oo. 

t&R J 

Hence, by the Prokhorov theorem, the family of measures {g^'^, t G M} is compact in T-L~^ for 
each e > 0. Therefore, the measures gf'"^ weakly converge to a limit measure g^^ as t — oo, 
i.e., for any e > 0, 



where g^'^ is a Gaussian measure in %. Therefore, gf — t g^ , t ^ oo, since gj^ = g^ ''^A'*. 
This implies convergence (Ill.Sp . ■ 

Remark The correlation matrix Q^{x, y) = (^0^(2; — 2/))ij=o,i) of the limit (Gaussian) measure 
g^ has a form 

q'c^{x-y) = ^(T++T_)f^(x-y), 

q'^ix-y) = -q'^ix-y) = ^{T+~T.)Vix-y), 

q]^{x-y) = \{T^+T^)5{x-y), 

where V{x) = —iFfT^^ [sgn(A;i)/u;(A;)]. In particular, the mean energy current density is for- 
mally 

v,-(o) = ^ vvm = -%=i^ f 45i^ * = — • (T. - r-, 0. . . . , 0). 



2{2nY J ^|A;|2 + m2 



Corollary 11.7 Let P : (0, ^) — )■ &e the projection of the space £^ on the space "H*. Let 
/If) be a Borel probability measure on such that the measure P/xq = I^-o — 9o > where g^ is 
constructed above. Then the measures fit weakly converge to a limit measure yUoo as t — )■ 00 on 
the space E'^ for each s < —d + 1/2. 
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11.4 Gibbs measure for coupled system 

Write (j) = (v?, vr) and C, = {q,p). For /3 > 0, we introduce the (normalized) Gibbs measures gjB 
on the space Formally, 



Z 



Definition 11.8 For P > 0, define the Gibbs measures gp{d(l)d^) in S^, s,a < —d/2, as 

g^{dcl>dO = ^e-^'-^^'^^^g^{d^) X g^idO- (11.10) 

Here (3 = 1/T is an inverse temperature, g^^dcj)) is defined in Definition \11.3l and g^ is the 
normalized Gibbs measure on x M'^, 

gM) = ^^(^) = l(\p\' + ^'I^H- (11-11) 

In Section n 1.61 we will prove the invariance of the Gibbs measures gjs w.r.t. the dynamics 

St. 

Lemma 11.9 The flow St is mixing w.r.t. g/s, i.e., for any functions Fi,F2 G L2{S^,gi3) we 
have 

j F,{StY)F2{Y)gp{dY) ^ j F,{Y)g^{dY) j F2{Y)g^{dY) as t ^ oo. 
Proof It is enough to check that for any Zi, Z2 G V, 

j {StYo, Z,){Yo, Z2) gp{dY,) ^0, t ^ 00. (11.12) 

Let Zi = {f,u,v) G v. By Eqns (13.11) . (13. 2p and by Theorem 13.11 we obtain 

\qt - {Wt(Po,(y)\ + \pt - {Wt^o, < C^em{t)\^o\+C2^/IJt)svip\{W^^o,'^po)\, 

and 

1(0^,/) - {Wt4>o,f.)\ < C^ernmo\ + C2^/IJJ)(snp\{Wr4>o,Vpo)\ + snp\{Wr(Po,a)\ 

These bounds can be proved similarly to Proposition 15. 1[ Hence, to prove (I11.12p it suffices to 
verify that 

j {(po,WU){Yo,Z2)gp{dYo)^Q as t ^ 00, (11.13) 

with X = G H'^ (see Definition 16.11) . Definitions 111.31 and lll.8[ equalities (lll.3p . 

formula (16.51) and the Lebesgue-Riemann theorem imply (I11.13p . ■ 



The results of next sections is a modification of the results of [2T], where the ergodic 
properties of a nonlinear system "heat bath-particle" were studied. 
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11.5 Effective Hamiltonian 

We rewrite the system fl3.ip - fl3.2p in new variables. Introduce an effective potential by 

V,ff{q) = ^{uj'\q\'-q-Kmq). (11.14) 

By condition Al', K=//(g) > 0. Define M'^-valued function h{x), x G M'^, 

h{x) = (A - m^)-^Vp{x), (11.15) 

wliere p is tfie coupled function, and put = {h, 0) E x R^. Then the fist equations in 
(12.1 p become 

4>t = CB(pt + qf'^P° = 'CB{4>t + qf ho), with Cb = ^ ^ ^2 (11.16) 

because £5/2.0 = (0, Vp). Define a M^-valued function ip = ipi^x), x G M*^, where 

^ = (^0,^1): ijj'^ = if + q ■ h, i/j^=7r. 

Then flll.l6p becomes ipt = C-B^Jt + ■ ho. Recall that £b is the generator of the group Wj. 
Hence, in new variables {ipt,^t) the system fl3.ip - fl3.2l) becomes 

t 

i)t = Wti^o + j Wt-sho ■ qsds, xER'', t G M, (11.17) 


t 

qt = -WV,ff{qt)- jv{t-s)q,ds + J'{t), (11.18) 



where J^{t) := (Vpo, W^tV'o), 

^V,ff{qt) = {uH - K^)qt, (11.19) 
the matrix is defined in (12.40 . its entries are 

J + m'' 

and r(t) stands for the M*^ x matrix with entries Tijlt), 

T,,{t):=-{V,po,WX) = i2n)-' I k,k,^^^^ \p{k)\' dk, z,j = l,...,d, (11.20) 

The equation flll.l8p is called the generalized or retarded Langevin equation with the random 
force T{t) and with the memory kernel T{t). 

Remark 11.10 (i) Eqn U1.20\) implies that r(0) = Km- Moreover, Tij{t) = —Dij{t), where 
Dijif) are the entries of the matrix D{t) (see liS. 3\) ). 

(a) By condition A2, po{x) = ho{x) = for \x\ > Rp. Hence, T{t) = for \t\ > 2Rp if m = 
(since in this case d>3 and odd), and \T(t)\ < C(l + |t|)~^/^ if m 7^ 0. 
(ill) J^{t) = F{t) - T{t)qo with F{t) from (EJ). 
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H Since j {F,{t)F,{s))g^pm = ( j (^(x) ® ^l^{y))g^{d^l^), (P^/V^po) (x) ® (VT/V.Po) (y)), 
then, by l\6.5\) . U1.3\) and ^11.2(]\) . we obtain 

This is an example of the fluctuation- dissipation theorem. 

Introduce an effective Hamiltonian H"^^ = + Hence, by fll.ip . 



Definition 11.11 (i) Define a map on S^: 

T:(</.,0->(^,0, iJ = ^ + q-ho. 



IS 



(a) Denote gj{dijjd^) := gp^T ^{dtpd^)). Then, gj{dtpd^) = g^{dilj) x g'^^-' ■' (d^) , where g^ i 
defined in Definition \11.3[ g'p^ is a Gaussian measure defined by g}" {d^ = ^e'^^A (O^^^. 

11.6 Invariance of Gibbs measures gj3 

Proposition 11.12 Let conditions A1-A3 hold. Then the Gibbs measures gp are invariant 
with respect to the dynamics, i.e. 

S;gp{uj):=gp{S;^uj) = gp{uj), for u e B{£'^) andteR. (11.21) 

Here B{£^) is the Borelian a -algebra of subsets in S^. 

Proposition 11.13 For every P > 0, the identity ^11.21\) holds. 

Proof. For simplicity, we omit indices a, s in notations £^ and "H^. The invariance f lll.2ip is 
equivalent to the identity: 

^ j F{StY)gp{dY) = 0, t G R, for every FiY) G Cb{8). (11.22) 

It suffices to prove (111.221) with t = only. Indeed, since St+s = St Sg, we have 

^ ^ FiSsStY)gpidY) = j-J^ F{S,+tY)gp{dY) = ^ ^ F{StSsY)gp{dY). (11.23) 



= for any fixed t G 

s=0 



Let i-J^F{S,Y)gp{dY) = 0. Then £ F{SsStY)g^{dY) 

s=0 

Hence, (I11.23P implies 

= 1/" F{S,StY)g^{dY) =^ [ F{StS,Y)gp{dY) =^ [ F{StY)g^{dY) 
as s=o at j£ s=o at 

and inT22^ follows. 
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Furthermore, it suffices to check flll.22p with t = and F{Y) = exp{i{Y, Z)) for every 
Z = (/o(x), A;, /) G D © D © M'^ © where D = C^{R'^). Then (0X221) with t = 

becomes formally (see f l2.3p ) 



Jtj/''''''''^9p{dY)\^^^ = jS^^'h{C{Y),Z)g^{dY) = 0, (11.24) 

where 

>C(¥?, vr, q,p) = (vr, (A - m'^)ip + q ■ Vp,p, -UJ^q + (Vp, y?)). (11.25) 
Now we prove flll.24p . Denote by / the integral 

I:= J^e^^''''h{CiY),Z)g^idY), 
and check that 7 = 0. Definition 111.111 implies g^ldY) = gJ{TdY). Hence, 

I F{Y)gp{dY) = I F{T~'Y)gJ{dY) = J gf{dO J F{^Ij - q ■ h,,i) g^^^. (11.26) 



v. 



Using flll.25p . (111.260 . and (Ill.lSp . we rewrite / in the form 



(11.27) 

+ik-p + il- [-u^q + (Vp, ipo-q- /i)] )^°(#o)^/^(#l)• 



w 



Integrals over Gaussian measures g'^{dipo) and gjj{dipi) can be represented as variational deriva- 
tives of their characteristic functional g^{fo) and 5'^(/i): 

,^(^'f)^{^,.)g^^{d^) = {±gl{f),.), ^ = 0,l, feC^{R'). 

Then 

6 



I (11.28) 
+zk-p + l-[- luS + Vp) - 2(Vp, q ■ h)\yg%f,)glUi) gf^- 
Using flll.4p we calculate 

{^am ■) = -ie-^«-^+™^)-^^'^) ((-A + m^)- V, ■) 
1 

Therefore, we reduce (I11.28P to the following integral 

I = C J e^e^-'J+'-rte-^^-i/^./o} (^,},.p_,i. VVeffiq) + (/o, /i)) e-^<'«) d^ 

= C^J e*(^-^+'-P)(fc . Vp - / ■ V,) [e-*«-<'''^«>-^^^"('''^')] dqdpgl{fo)gl{h), 

by (111.140 . (Ill.lSp and (111.190 . Partial integration in q and in p leads to 

I = C2 [ e^(^-«+'-P)(-A; ■ (il) + / ■ (ik))e~'''-'^f'^M-P"A"ii^p) dqdp = 0. ■ 
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